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The cosmic distance duality relation is a milestone of cosmology involving the luminosity and an¬ 
gular diameter distances. Any departure of the relation points to new physics or systematic errors 
in the observations, therefore tests of the relation are extremely important to build a consistent 
cosmological framework. Here, two new tests are proposed based on galaxy clusters observations 
(angular diameter distance and gas mass fraction) and H{z) measurements. By applying Gaus¬ 
sian Processes, a non-parametric method, we are able to derive constraints on departures of the 
relation where no evidence of deviation is found in both methods, reinforcing the cosmological and 
astrophysical hypotheses adopted so far. 

PACS numbers: 


I. INTRODUCTION 


One of the key relations in cosmology is the cosmic 
distance duality relation (CDDR), which expresses a con¬ 
nection between two observables: the luminosity distance 
(11 and the angular diameter distance c?a through 


dh ^ ^ .. 

{l + zfdA ’ 


( 1 ) 


where z is the redshift of the source. 

Although it appears to be a fortuitous relation in an 
FLRW (Friedmann - Lemaitre - Robertson - Walker) 
universe in many textbooks, it was actually derived by 
Etherington long ago [I|, where its validity implies pho¬ 
ton number is conserved, gravity is described by a metric 
theory and light propagates on unique null geodescis Q . 
Therefore, any violation of this relation is a clear indica¬ 
tion of new physics. 

Approaches to test the CDDR in the recent literature 
assume a cosmological model suggested by a set of ob¬ 
servations (usually the ACDM concordance model) and 
check the validity of the CDDR in the context of some as- 
trophwical effect. Examples of these methods are given 
bylli. Other tests have been proposed with the goal 
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of testing effectively such relation, involving from lumi¬ 
nosity distances of type la supernovae, angular diameter 
distances of galaxy clusters, gas mass fractions of galaxy 
clusters, the blackness of the cosmic microwave back¬ 
ground and gamma-ray bursts [e.g. 0,0,0,0,0,0,[10,[ll| . 
Generally, departures from the relation are written as 
r]{z) = (1 -b z)~'^dL/ cLa, where for 77 is adopted simple 
parametric forms as r]{z) = 1 - 1 - rjoz or rj^z) = 1 - 1 - 770 

An interesting test for the distance duality relation 
was recently proposed by [l 0 ], where one can write 
V = f Sze/ f x-ray, fszE Standing for the gas mass frac¬ 
tion of a galaxy cluster measured through the Sunyaev- 
Zel’dovich effect [T 0 | and f x-ray for the gas mass frac¬ 
tion obtained calculating the gas mass observed in X-rays 
divided by the total mass assuming hydrostatic equilib¬ 
rium. The method has some advantages comparing to 
other approaches: as the same object is used to deter¬ 
mine 77, there are no problems of comparing distances 
in completely different regions of the sky where inho¬ 
mogeneities can be important, as well as it is easier to 
grasp systematic errors. Considering the two standard 
parametrizations for 77 described above, no violation for 
the relation was the derived for a subset well described 
by hydrostatic equilibrium. 

In this paper, two new tests are proposed based 
on the distance to galaxy clusters from the ESZ/X- 
ray technique. X-ray gas mass fraction and H(z) mea¬ 
surements. The use of H{z) measurements allows a 
model-independent test, overcoming difficulties when us¬ 
ing galaxy cluster and type la supernovae (SNe la) due 
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to a dependence on a cosmolo gica l model and the Hubble 
constant when using SNe la [ij. Also, another advan¬ 
tage of our tests is that they rely on data coming from 
the observation of the same galaxy clusters, which allow 
us to search for new physics as well as to check the con¬ 
sistency of the hypotheses adopted in such observations. 
As we shall see, the results are derived using Gaussian 
Processes, a non-parametric method which does not re¬ 
strict a functional form for rj[z), where both methods 
fully agree with the standard value 77 = 1 in all redshift 
range covered by data. The paper is organized as fol¬ 
lows: in Sec. El we describe the sample used in the work. 
Sections imi and m details the methods and the non- 
parametric approach adopted, respectively. The results 
are in Sec. |Vl We close the paper in Sec. I^with the 
conclusions. 


II. SAMPLES 


In order derive constraints for ry free from parametriza- 
tions we use: 

(a) 38 gas mass fraction of galaxy clusters from [Tsj 
spanning redshifts from 0.14 to 0.89, derived from Chan¬ 
dra X-ray data. These authors used the non-isothermal 
double /3-model for gas distribution. The 3D tempera¬ 
ture profile is modeled assuming that the ICM is in hy¬ 
drostatic equilibrium with a NFW dark matter density 
distribution (T^ . 

(b) 38 angular diameter distance from galaxy clusters 
obtained from their Sunyaev-Zeldovich and X-ray obser¬ 
vations (the so-called SZE/X-ray technique). The sample 
is in redshift range is 0.14 < 2 : < 0.89 and was compiled 
by [13 where the cluster plasma and dark matter distri¬ 
butions were analyzed assuming a non-isothermal spher¬ 
ical double /3 model. T he g alaxy clusters of this sample 
are exactly the same of [l^ . 

(c) 28 Hubble parameter versus redshift data points, 
H{z), from cosmic chronometers, BAOs and a local mea¬ 
surement in redshifts up to 2 = 1.75 [H, [l^ HO, Hlj, 1^ . 

For a compilation of H{z) measurements see e.g. 

13. 


III. METHODS 


Method I 


The gas mass fraction of galaxy clusters is given by 

M 


fgas — 


Mga, 

Mtot 


oc dhdA^^'^ , 


( 2 ) 


where Mgas is the gas mass and Mtot is the total mass 
(including dark matter) of the galaxy cluster. As it is 
known, along with the assumption of a constant fgas with 
redshift it is possible to test different cosmological sce¬ 
narios. However, in a recent paper, [0| showed that the 
gas mass fraction is dependent on ry value. These authors 
deduced a more general expression for the gas mass frac¬ 
tion given by 


fgas{z) = N 


dLdT 


( 3 ) 


By using a general version of equation (1) (?y ^ 1) was 
obtained: 


dA{z) = 


d\ 




( 4 ) 


where N is the normalization factor that carries all the 
information about the matter content in the cluster and 
the quantities with the subscript * correspond to those 
used in the observations (concordance model where rj = 
1). In this way, it is possible to define 


Vobsiz) = N 


d 


*3/2 

A 


d^/^fgasiz) 


( 5 ) 


In our analysis, dA is obtained directly from H{z) mea¬ 
surements, as it will be explained later. It is important to 
stress that the test for ry proposed here is independent of 
the assumption of a constant fgas with redshift since we 
are interested in dA from each cluster separately and not 
from its distribution over 2 . The sources of uncertainty 
in the measurement of fgas to the galaxy cluster samples 
are: i) the statistical errors concerning to the X-ray and 
SZE observations comprehend known effects: SZE point 
sources ±4%, kinetic SZ ± 8 %, ±20% and ±10% for clus¬ 
ter asphericity to fgas from X-ray and SZE observations, 
respectively, ii) the systematic errors for the galaxy clus¬ 
ters are: X-ray absolute flux calibration ± 6 %, X-ray tem¬ 
perature calibration ±7.5%, SZE calibration ± 8 % and a 
one-sided systematic uncertainty of — 10 % to the total 
masses which accounts for the assumed hydrostatic equi¬ 
librium (see subsection 5 in (l^l)- In our analysis the 
systematic errors will be add in quadrature to statistical 
errors [T7I |. 


While the item (c) in the earlier section is insensitive 
with respect the cosmic duality relation validity, the as¬ 
tronomical observations in items (a) e (b) are dependent 
through different ways. In our analysis we use two meth¬ 
ods to obtain constraints on ry, namely: 


Method II 

As it is largely known, the X-ray emission of the in¬ 
tracluster galaxy medium can be combined with the 
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Sunyaev-Zeldovich effect to estimate the angular diame¬ 
ter distance (ADD). The main advantage of this method 
is being fully independent of any local calibrator. In the 
context of this phenomenon, the angular diameter dis¬ 
tance is such that 


dA\ATo)^A eHO 1 ^ (Aro)^Aego 1 
dL'^SxoTeo'' 00°" {I + z)^SxoTeO^ ’ 


( 6 ) 


where Sxo is the central X-ray surface brightness, Teo 
is the central temperature of the intra-cluster medium, 
AeHo is the central X-ray cooling function of the intra¬ 
cluster medium, ATq is the central decrement tempera¬ 
ture, and 6c refers to a characteristic scale of the cluster 
along the line of sight (l.o.s.), whose exact meaning de¬ 
pends on the assumptions adopted to describe the galaxy 
cluster morphology. 

In Q it was shown that this technique is heavily depen¬ 
dent on the cosmic distance duality relation. If one as¬ 
sumes a deformed expression such as equation (I) (77 7 ^ I) 
it is possible to obtain that 

dr^‘^{z)=dA{z)v{z?. (7) 


This quantity is reduced to the standard angular diam¬ 
eter distance only when the CDDR relation is strictly 
valid (77 = 1). In this way, it is possible to define 


Vobsiz) 


/^data(^)\ 1/2 

V dAiz) ) 


( 8 ) 


In our analysis, dA{z) is obtained directly from H{z) 
measurements. The sources of uncertainty in the mea¬ 
surement of to the galaxy cluster samples are: i) 

statistical contributions: SZE point sources ±8%, X-ray 
background ±2%, Galactic Nh < ±1%, ±15% for clus¬ 
ter asphericity, ±8% kinetic SZ and for CMB anisotropy 
< ±2%, ii) estimates for systematic effects are as follow: 
SZ calibration ±8%, X-ray flux calibration ±5%, radio 
halos -1-3% and X-ray temperatute calibration ±7.5% (see 
table 3 in [13 )■ 


Distances from H(z) measurements 


The H(z) measurements were derived from the follow¬ 
ing physical observables: cosmic chronometers [3, [l^, 
I21I . [23, baryon acoustic oscillations obtained from the 
clustering of galaxies [l^ [13 and from the Wide Field 
Camera 3 (WFC3) on the Hubble Space Telescope (HST) 


The method of cosmic chronometers [23 is based on 
the relative age of a galaxy population at different red- 
shifts. As H{z) = for spectroscopic redshifts 

with interval A 2 ;, a measured age difference At provides 
a direct inference of H{z). Note that this measurement 
is independent of absolute ages as well as the metric con¬ 
sidered. The error budget considered is a combination 
of statistical and systematic errors, whereas the former 


is due to the calibration relation to infer the ages and 
the latter is due to metallicity, stellar population synthe¬ 
sis model and star formation history. Other sources of 
systematic errors are not relevant given the size of the 
other errors, as progenitor bias, initial mass function and 
a-enhancement. 

The values of H{z) derived by [13 come from the 
baryon acoustic oscillations (BAOs), while the ones ob¬ 
tained by [13 come from combining the BAOs with the 
Alcock-Paczinsky test [13- The BAOs come from the 
sound waves propagating in the photon-baryon fluid in 
the early Universe. After the decoupling, these modes 
were frozen and a preferred scale was imprinted on the 
distribution of matter, which can be detected analyzing 
the distribution of objects as galaxies. This scale is the 
sound horizon at the baryon drag epoch, so it can be 
determined from observations of the cosmic microwave 
background radiation. By measuring this scale, which is 
a standard ruler, at transversal and radial directions, one 
is able to derive the angular diameter distance and the 
Hubble parameter, respectively. The Alcock-Paczinsky 
test relies on comparing the radial and transversal sizes of 
objects assumed to be isotropic. Equating their sizes, one 
is able to derive constraints on a combination of angular 
diameter distance and Hubble parameter independently 
of a standard ruler. Both tests are weakly dependent on 
the adopted cosmological model, therefore it can be used 
in cosmological tests. [l3 estimated the impact from 
systematic errors in the measurement of the correlation 
function as well as the BAG peak location, where simula¬ 
tions were extensively used. Possible sources of system- 
atics include the choice of radial selection function, finite 
volume effects and the accuracy of the statistical error 
model. In this work we considered the statistical and 
systematic errors obtained in [13 ■ [13 discarded pos¬ 
sible sources of systematic errors as much smaller than 
statistical errors, as modelling redshift-space distortions, 
varying fitting range for the BAO peak and implemen¬ 
tation of a quasi-linear model. Also, comparison with 
constraints coming from simulations found no evidence 
of systematic errors. 

We follow the methodology firstly presented by [13 
where one transforms H{z) measurements into cosmo¬ 
logical distance estimates by solving numerically the co¬ 
moving distance integral for non-uniformly spaced data, 
i.e.. 


Dc 


r dz' 

Jo 'W) 


c 

2 


N 


1 


1 

H{z^)\ ’ 


( 9 ) 

Since the error on 2 ; measurements is negligible, we only 
take into account the uncertainty on the values of H{z) 
(the error from the method itself is also completely negli¬ 
gible compared to the errors of H{z)). As one may check, 
by using standard error propagation techniques, the error 
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associated to the bin is given by 


1/2 


Si = - ^i) 




"i+i 


Hf 


( 10 ) 


so that the error of the integral m in the interval z = 0 
- Zn is Si- Thus, we use 28 independent H(z) 

measurements [3, [3, [13, EH [H, EH EH in the redshift 
range 0 < z < 1.75. 

As one may see, the comoving distance obtained from 
this technique is completely independent of the cosmic 
distance duality relation. We restrict our analysis to the 
flat case, in this way the ADD in methods (I) and (II) is 
cIa = (I + z)~^Dc- Moreover, since the galaxy clusters 
and H{z) observations are performed at different z, we 
calculate Dc at each galaxy cluster redshift from a poly¬ 
nomial fit of the Dc from H{z) points shown in Fig. (1). 
In Fig. (2) we show the Tjobs values for each method by 
using equations (5) and ( 8 ). For method (I) we marginal¬ 
ize over the parameter N. It is important to stress that 
all systematical errors are being considered. 


IV. GAUSSIAN PROCESSES 

Gaussian Processes (GPs) is a non-parametric method 
which allows one to reconstruct a function without spec¬ 
ifying a form for it. GPs generalize the notion of a 
gaussian random variable for a random function, being 
characterized by a mean and a covariance function. Our 
knowlegde about the properties of the function are given 
in the covariance function, where we adopt the standard 
squared exponential covariance function k{z,z): 



Figure 1: The black squares points are the comoving distance 
obtained from H(z) measurements. The red line solid corre¬ 
spond to polynomial fit and the dashed lines the Icr error. 



k{z, z) = Of exp 


2/2 


( 11 ) 


where ct/ and I are hyperparameters which control how 
the function changes in the y and x axes, respectively, 
and the covariance function relates to points in input 
space z and z. By considering that the function we want 
to reconstruct is a given realization of a GP, we have to 
calculate the probability of observing a new set of points 
at z* with mean values f* and covariance cov{t*) given 
we measured y at positions z with covariance matrix G, 
where the mean and covariance are given by 


f* =K{z*,z)[K{z,z)+C]-^y, (12) 


cov{r) = K{z*,z*) - K{z*,z)[K{z, z) -b C]-^K{z, z*). 

(13) 

In the above expressions Ar(zi,Z 2 ) is a matrix connect¬ 
ing all possible values of Z\ and Z 2 . The hyperparameters 
must be margmalized over, where we follow the approach 
described in |3l| . GPs have been adopted in many cosmo¬ 
logical applications [e.g. EH Eli EH EE Eli EH EE El) nil ■ 


Figure 2: y{z) values. The filled red circles and the open 
blue circles are the results from Method I and II, respectively, 
considering all systematic errors. 


V. RESULTS 


In order to derive constraints to ri{z) two steps are 
needed. Firstly, we derive 77 for each cluster and its 
error for Methods I and IF From these values we are 
able to reconstruct r]{z) by applying a squared exponen¬ 
tial covariance function and imposing the condition that 
77 ( 0 ) = 1.0. Figure |3] and Fig. |4]show the constraints on 
r]{z) for methods I and II, respectively. The shaded con¬ 
tours represent 68 % and 95% confidence regions and all 
known sources of statistical and systematic errors were 
taken into account in our analyses. For comparison it is 
also plotted the results of a linear fit ri{z) = 1 + tjqz with 
the 95% confidence levels as red dashed lines. 

As one can see, both methods are compatible within 
the 68% conhdence level with 77 = 1 , i.e. with the validity 
of the CDDR. Apart from not showing any evidence for 
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Figure 3: Constraints to rj{z) from Method I. The black 
line corresponds to the mean reconstructed function and the 
shaded contours 68 % and 95% confidence levels. All known 
sources of systematic errors were included in the analysis and 
the standard value r] — 1 is within the 68 % for the whole red- 
shift region. The red dashed lines represent the mean value 
and 95% conhdence levels of a linear fit 77 ( 2 ) = 1-1- rjoz. 



z 


Figure 4: Constraints to r]{z) from Method IF As in Fig. 
[ 3 ] the black line shows the mean reconstructed function and 
the shaded contours represent 68 % and 95% confidence lev¬ 
els, where the red dashed lines depict the best ht and 95% 
confidence levels of the linear ht. Again all known sources 
of systematic errors were considered and 77 = 1 is within the 
68 % conhdence level. Method II provided smaller errors than 
Method I, although both methods show no deviation of the 
cosmic distance duality relation. 

new physics, our results also point to an internal consis¬ 
tency for deriving constraints from galaxy clusters with 
the gass mass fraction and with the angular diameter 


distance. 

In order to see whether the results were consistent 
given the adopted method, we redid our analysis as¬ 
suming different covariance functions for the GP regres¬ 
sion. We considered some instances of the Matern family, 
see [13 and [13 for definitions and discussion. The re¬ 
sults were practically indistinguishable from the standard 
squared exponential covariance function, which reinforces 
the robustness of the results. 


VI. CONCLUSIONS 

In this work we have showed that observations of 
galaxy clusters jointly with cosmic expansion rate can 
be used to derive constraints to the distance duality re¬ 
lation (CDDR), = ?7 = 1 . From galaxy clusters 

observations we used, separately, a sample of 38 angular 
diameter distances obtained via their gas mass fraction 
in X-ray band [l5l| and another one of 38 angular diam¬ 
eter distances obtained via their Sunyaev-Zeldovich ef¬ 
fect and X-ray observations [I3- The galaxy clusters are 
the same in the samples and they are in redshift range 
0.14 < 2 < 0.89. As stressed in this paper, both obser¬ 
vations are heavily dependent on the CDDR validity. In 
order to obtain CDDR independent distances and per¬ 
form our analyses we have used 28 recent H (z) measure¬ 
ments. The H{z) measurements are transformed into 
cosmological distance estimates by solving numerically 
the comoving distance integral for non-uniformly spaced 
data (see Fig. (1)). Unlike what happens in most papers 
in the literature no parametrization was used to 77 , where 
we used Gaussian Processes, a non-parametric method 
to reconstruct the 77 function. As a general result we ob¬ 
tained that validity of the relation was ensured within 
1(7 for the redshift range 0.14 < z < 0.89 including all 
systematic errors to both galaxy clusters observations. 
Since some Sunyaev-Zeldovich surveys are underway, or 
planned to start in the near future, more and larger data 
sets with smaller statistical and systematic uncertainties 
will be available and the method proposed here may im¬ 
prove the constraints on fundamental physics as well as 
the astrophysics of galaxy clusters. 
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